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We propose a simple method for incorporating correlations into the impact parameter space 
description of multiple (semi-)hard partonic collisions in high energy hadron-hadron scattering. 
The perturbative QCD input is the standard factorization theorem for inclusive dijet production 
with a lower cutoff on transverse momentum. The width of the transverse distribution of hard 
partons is fixed by parameterizations of the two-gluon form factor. We then reconstruct the hard 
contribution to the total inelastic profile function and obtain corrections due to correlations to the 
more commonly used eikonal description. Estimates of the size of double correlation corrections are 
based on the rate of double collisions measured at the Tevatron. We find that, if typical values for 
the lower transverse momentum cutoff are used in the calculation of the inclusive hard dijet cross 
section, then the correlation corrections are necessary for maintaining consistency with expectations 
for the total inelastic proton-proton cross section at LHC energies. 
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I. INTRODUCTION 

Models of multiple partonic collisions in high energy 
hadron-hadron scattering are important for simulations 
of complex events in upcoming experiments like those 
at the LHC or in high energy cosmic ray air showers. 
Furthermore, measurements of the rate of multiple par- 
tonic collisions can be used to test current models of pro- 
ton structure. Already, measurements at accessible en- 
ergies [H-d] yield a much smaller effective cross section 
CTcff than what is naively expected if partons are homo- 
geneously distributed over the transverse area of the pro- 
ton. (Figure [T] shows a schematic depiction of a double 
partonic collision.) The definition of Coff is 
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where a2 is the inclusive cross section for a single par- 
tonic collision (resulting in a dijet), tT4 is the inclusive 
cross section for a double collision, and m is a symmetry 
factor that depends on whether the partons are identi- 
cal. New measurements of multiple collisions are cur- 
rently being proposed for the LHC [6]. Hence, novel new 
phenomena involving multiple hard partonic collisions, 
which will enhance understanding of proton structure, 
can be expected in the next generation of experiments at 
the high energy frontier. 

However, multiple interactions involve a complex inter- 
play of soft, hard, and semihard physics, so a complete 
description using purely perturbative techniques is not 
possible. In simulations of complex hadronic final states, 
methods are needed for combining hard and soft colli- 
sions in a consistent way. (For an overview of current 



approaches see, e.g. 0,(1], and references therein.) The 
hard contribution is calculated using the well-known per- 
turbative QCD (pQCD) leading twist factorization for- 
mula for the inclusive dijet production cross section, in- 
volving a convolution of the standard parton distribution 
functions (PDFs) with a partonic cross section. An im- 
mediate complication is that, in order for perturbation 
theory to be applicable, the relative transverse momen- 
tum of the produced jets must be larger than some mini- 
mum cutoff scale p^. The cutoff should be chosen small so 
as to maximize the range of the perturbative expression, 
but still large enough for perturbative methods to be rea- 
sonable. For describing events with transverse momen- 
tum less than p^, nonperturbative methods are needed. 
A prescription for matching the high and low transverse 
momentum behavior is necessary for a complete descrip- 
tion over the full range of transverse momentum. 

The inclusive pQCD hard dijet cross section is numer- 
ically very sensitive to the precise choice oip^. This has 
a tendency to lead to substantial variations between dif- 
ferent model predictions of the minijet cross section at 
high energies. The question of what values of pi are ap- 
propriate continues to be discussed in current research 
on the development of models and simulations (see, for 
example, recent discussions in [st-flTj). 

An additional complication is that a description is 
needed for the distribution of hard partons in impact 
parameter space. In the past, it has usually been mod- 
eled or assumed to be equal to the electromagnetic form 
factor of the nucleon. Also, in a number of models, soft 
and hard partonic interactions are incorporated into a 
single eikonal picture. In this way, both soft and hard in- 
teractions are included in a way that respects s-channel 
unitarity. Data for the total and inelastic pp cross sec- 



FIG. 1: Schematic depiction of a double hard coUision. A disconnected pair of partons from each proton colhde to produce a 
pair of high transverse momentum dijets. 



tions can then be used to fit parameters such as pi and 
the width of the distribution of hard partons in impact 
parameter space. However, different choices for these pa- 
rameters can provide equally good fits to the total cross 
section at accessible energies while leading to very dif- 
ferent extrapolations at high energies (see, for example, 
Ref. [13 )■ Therefore, it is important to make use of any 
experimental or observational information that can nar- 
row the range of allowed parameters and falsify some of 
the models currently in use. 

It is nowadays possible to use pQCD to obtain direct 
experimental information about the transverse spread of 
hard partons in the proton via parametrizations of the 
generalized parton distributions (GPDs). (See Ref. [ll] 
and references therein for a review of the phenomenology 
of GPDs.) The gluon GPD of the proton, for example, 
can be extracted from measurements of the t dependence 
in deep inelastic production of light vector mesons or 
photoproduction of heavy vector mesons [3, [3 ■ Then, 
since the GPD is a universal object [11], it can be reused 
in other processes. In particular, it can be used in the 
description of the impact parameter dependence of hard 
collisions in hadron-hadron scattering. As such, we adopt 
the point of view in this paper that the impact parameter 
description of a hard collision is not an adjustable model 
parameter, but rather is fixed by fits to the gluon GPD. 

An accurate description of the impact parameter de- 
pendence is important because it allows for a measure 
of how close the hard interaction is to the unitarity (or 
"saturation") limit. Furthermore, the proximity to the 
unitarity limit is related to the number of high transverse 
momentum jets that are produced, and hence to the use 
of central! ty as a dijet trigger .17] . However, as was il- 
lustrated in Ref. [l^l, using GPDs to describe the im- 
pact parameter dependence of hard collisions within the 
approximation where partons are not correlated in the 
transverse plane (which is implemented in the simplest 
eikonal description) leads to inconsistencies with general 
expectations for the extrapolation of the total inelastic 
profile function to high energies, unless a large value is 
used for the transverse momentum cutoff p^. Specifically, 



the contribution to the total inelastic profile function 
from the production of hard dijet pairs becomes larger 
than the total inelastic profile function itself. Although 
the inclusive dijet cross section is unitarized in the ba- 
sic eikonal description, it nevertheless grows too rapidly 
with energy. As we will argue, this problem is most likely 
a symptom of the common assumption that partons are 
uncorrelated. 

The study of correlations in multiple hard collisions 
is already an active area of research (l9l - [23 . Imple- 
mentations of Dokshitzer-Gribov-Lipatov-Alterelli-Parisi 
(DGLAP) evolution in multiparton distribution func- 
tions suggest that correlations are indeed significant [2J . 
Correlations may also be induced by evolution in Bjorken 
X psj in the very high energy limit. However, numerical 
estimates in Ref. [18,] suggest that hard unitarity (satu- 
ration) effects contribute to only a small fraction of the 
total inclusive dijet cross section, even at LHC energies. 
Moreover, the inconsistencies encountered in Ref. [18] oc- 
cur even at large impact parameters, > 1.0 fm, and values 
of that are not particularly small, e.g. p^ « 2.5 GeV. 
In the instanton liquid model (see Ref. !24] and references 
therein) one can expect a strong correlation between the 
quark and gluon fields in the constituent quarks which 
will result in correlations at all impact parameters. It 
may be somewhat diluted at large b due to the contri- 
bution from the 4gg component. On the other hand for 
partons at large enough distances from the center it is 
likely that another parton should be present at a distance 
comparable to the confinement scale, leading to correla- 
tions of qq pairs at relative distances < 0.5 fm that are 
rather small which could easily compensate for the dilu- 
tion effect. Therefore, it is likely that nonperturbative 
correlations, unrelated to the overall impact parameter, 
play a role in determining the size of the dijet contri- 
bution to the total inelastic cross section. Heuristically, 
one can imagine scattering at large impact parameters 
as the scattering of pion clouds. The nonperturbative 
dynamics responsible for binding the qq pairs introduces 
potentially large correlations. 

CoUinear constituents of the incoming protons are ex- 
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pected, on average, to have a spacetime separation of or- 
der ~ 1/Aqcd. Because of the breakdown of asymptotic 
freedom at these scales, it should be anticipated that the 
constituents of the nucleon are subject to strong nonper- 
turbative correlations. In other words, one should not 
expect the multiparton distribution functions to be sim- 
ply related to products of the single parton distribution 
function at any scale. 

As a specific illustration, let us consider the double 
parton event in Fig. [T] If /iv(2;i,a;2) is the distribution 
function for parton pairs with momentum fractions Xi 
and X2 inside proton iV, and fN{xi) and /at (0:2) are the 
standard PDFs, then an often made ansatz is 

?? 

fN{xi,X2) fN{xi) fN{x2)- (2) 

In some treatments of correlations this relation follows 
after an integration over impact parameters, even if the 
factorization ansatz is broken in impact parameter space. 
The question marks in Eq. ^ indicate that we will ques- 
tion the validity of this assumption. Given the strong 
binding between the constituent partons in the proton, 
the approximation in Eq. ^ is probably very rough at 
any scale. We will argue that deviations can have an 
important effect on extrapolations to high energies. 

For a large cutoff , the inclusive dijet cross section 
gives a small fraction of the total inelastic cross section, 
and there is no conflict with s-channel unitarity. What 
is needed at smaller pi is a method for organizing corre- 
lation corrections which does not rely on the commonly 
used assumption that correlations can be neglected. 

To summarize, the aim of this paper is to set up a 
method for organizing corrections to the uncorrelated ap- 
proach at small p^. We will use this to simultaneously 
incorporate the following information into a description 
of multiple hard partonic collision: 

• The impact parameter distribution of hard partons 
obtained from measurements of the gluon GPD. 
The width of the distribution of hard partons is 
not a fixable parameter in our approach. 

• An estimate of the size of double parton correla- 
tions obtained from measurements a^s. In this 
paper, we will assume that these correlations are 
roughly independent of impact parameter. 

From this information, we will reconstruct the hard dijet 
contribution to the total inelastic profile function. Com- 
pared with the uncorrelated expression, the result ob- 
tained with correlations will be shown to exhibit greater 
consistency with common high energy extrapolations of 
the total inelastic cross section, even with a relatively 
small and fixed value for p^ in the perturbative inclusive 
dijet calculation. Thus, including correlations in this way 
may provide a natural resolution to the consistency prob- 
lems encountered in Ref. [l8| at large impact parameters 
and fixed p^. 

The paper is organized as follows: In Sec.|ll]we review 
the basic setup for discussing high energy collisions in im- 
pact parameter space, and in Sec. lIIII we review the steps 




FIG. 2: Momentum labels for pp scattering. 

outlined in Ref. for dealing with multiple hard colli- 
sions in terms of inclusive dijet cross sections, applying 
the result to the special case of uncorrelated scattering in 
Sec. HVl In Sec. |V] these steps are extended to allow for 
impact parameter independent correlations. In Sec. IVII 
we use estimates of correlations based on Tevatron mea- 
surements of (Jeff to calculate the hard contribution to the 
inelastic profile function at LHC energies, and we com- 
pare with standard extrapolations of the total inelastic 
profile function. We speculate on prospects for including 
impact parameter dependence of correlations in Sec. lVIll 
We close in Sec. IVIIII with a conclusion and a discussion 
of the main results. 



II. IMPACT PARAMETER REPRESENTATION 

A. S-Channel Unitarity 

We work in impact parameter space, defining the pro- 
file function, 

r(.,5).^/d^qe--A(M), (3) 

Here A{s, t) is the amplitude for elastic pp scattering (see 
Fig. [2]), and s and t are the usual Mandelstam variables, 
s = {Pi + ^2)^ and t = {pi —psY- We work in the 
high energy limit, s » —t where one may approximate 
t fa — (see e.g. Ref. [2^ for a review of kinematics 
in the high energy limit). In the two-dimensional Fourier 
transform to coordinate space, b is the impact parameter. 

Unitarity and analyticity allow the total, elastic and 
inelastic cross sections to be calculated in terms of the 
profile function via the familiar relations: 

cTtot{s) = 2 J d%Rcris,b), (4) 

aei{s) = Jd^h \T{s,b)\\ (5) 

a„,,i{s) ^ J d^h (2ReTis,b) - \T{s,b)\^y (6) 

We refer to the integrand of Eq. ([6]) as the inelastic profile 
function, 

r->(s,6)^(2Rer(s,6)-|r(s,6)|^). (7) 
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In the very high energy hmit, it is appropriate to neglect 
the imaginary part of the amphtude. Then unitarity re- 
quires 



^incl 



(s,6),r(s,fe) < 1. 



(8) 



The profile function for the total proton-proton cross sec- 
tion and its extrapolation to high energies has been stud- 
ied extensively over the last few decades. Our main con- 
cern in this paper is whether common extrapolations of 
Eq. ([7]) are consistent with extrapolations of the total 
hard contribution. 



B. Inclusive Dijet Cross Section 

Hard collisions are described by the leading twist 
pQCD expression for the inclusive dijet cross section, 



i,j.k.l 



da. 



ij-^kl 



fi/pi{xi;pt) fj/p2{x2;pt)0{pt - Pt) ■ (9) 



The collision is between parton types i and j inside pro- 
tons pi and p2 respectively, and the partonic hard scatter- 
ing cross section doij^u I dpi is calculated at tree level. 
We have explicitly included a K factor and a symme- 
try factor 1/(1 -f Ski). The fi/p,{xi;pt) and fj/p^{x2;pt) 
are the usual integrated parton distribution functions for 
partons with longitudinal momentum fractions xi and 
X2, evaluated at a hard scale equal to the dijet trans- 
verse momentum. The lower bound pi is in principle 
arbitrary, but it should be chosen large enough for it 
to be a reasonable hard scale. It is not clear exactly 
what is the minimum p^ that can be used, but Relativis- 
tic Heavy Ion Collider (RHIC) data for pion production 
suggest that perturbation theory is still reliable for pi- 
ons with > 1 GeV [l^ for forward production where 
background from soft physics is small. (Note here that 
for such kinematics pt of the progenitor quark is close to 
Pt of the pion.) 

Information about the impact parameter distribution 
of hard partons in the proton is obtained from the gluon 
GPD which is parametrized in experimental measure- 
ments of the t dependence in exclusive heavy vector me- 
son photoproduction and light vector meson electropro- 
duction at small x. For the t dependence of the differen- 
tial exclusive vector meson production cross section, we 
use the dipole parametrization. 



Fg{x,t\y) 



1 



(10) 



The basic Feynman diagram contributing to this reaction 
involves the exchange of two gluons in the t channel, so 
we refer to it as the two-gluon form factor. The gluon 
GPD g{x,t;fi) is then 



g{x, t; fi) = g{x; y)Fg{x, t; /i), 



where g{x; ^) is the standard integrated gluon distribu- 
tion function evaluated at hard scale fi. The two-gluon 
form factor obeys the condition. 



Fg{x,t^O;fi) = 1 



(12) 



so Eq. ([TT|) reduces to the standard PDF in the limit of 
t — 0. The parameter mg(a;,/j,) in Eq. (|10p determines 
the width of the peak around t = 0. Following Ref. jT7| . 
we allow it to have x and fj, dependence to account for 
evolution in the hard scale fi and diffusion at small x. 
The Fourier transform of the two-gluon form factor into 
the transverse plain is 



J^g{x,b;fi) = J d^ At Fg{x,t; ^)e~ 



(13) 

Again, the approximation t « — At is justified so long as 
s >> —t. Using the dipole form in Eq. ([T0|) one finds 
explicitly 

Tgix,b;^,) = "''^''^f'\ ,imgix;fi)b). (14) 

In this paper, Kn for integer n denotes a modified Bessel 
function of the second kind. The overlap function is de- 
fined as 

P2{b,xi,X2;y) = 

j d^h' Fg{xl^h'\■^l)Fg{x2^h-h'\■^l). (15) 

Using Eq. yields 



P2{sMp'i) 



"^g(a;;Pt) / mg{x;p1)b 



127r 



K3{mg{x;p1)b). 

(16) 

Here we have made the usual approximation, xi ~ X2 ~ 
X = 2pj/i/i. A more detailed treatment should take into 
account the separate integrations over xi and X2 — there 
is not, in general, a one-to-one mapping between values 
of Pt and xi{2). For now, we mention that direct numeri- 
cal calculations verify that this approximation introduces 
less than 10% error in the essential region of integration 
for the cross section. Note that P2{s,b;p'j:) is normalized 
to unity. 



J d'h P2{s,b;p1) 



1. 



(17) 



Combining Eqs. (|10j|17p with Eq. © allows the inclusive 
dijet cross section to be written in the form 



d^hx2{s,b;p't) 



where 



(18) 



(19) 



(11) 



X2(s,6;pt) = c7'^^cD{s-:Pt)P2is,b;pt). 

We will refer to X2{s, b]Pi) as the impact parameter de- 
pendent inclusive cross section j4j]. Using the GPD 
to write the inclusive dijet cross section in the form of 
Eq. ((T5)) enables one to analyze the contribution from 
different regions of impact parameter space. 
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III. MULTIPLE HARD PARTONIC COLLISIONS 

By taking into account multiple hard scattering events, 
it is in principle possible to reconstruct the total hard 
scattering contribution to the total inelastic profile func- 
tion using probabilistic arguments [13, HI] . In this sec- 
tion, we briefly review the steps for constructing the di- 
jet contribution to the total inelastic profile function in 
terms of a series involving the inclusive n-dijet cross sec- 
tions. 

We start by defining X2n(s, to be the analogue of 
X2(s, b;Pf) for the case of rt-dijet production {n hard col- 
lisions). Namely, integrating over all impact parameters 
yields 

J d^hx2nis,b;p^,)^aZ{s;p'i), (20) 

where o'Z{s]p'[) is the inclusive cross section for produc- 
ing n-dijet pairs. 

Next, X2n(s,6;Pt) is defined to be the exclusive ana- 
logue of X2n(s,&;Pt). It describes the production of ex- 
actly n hard collisions, differential in b. Integrating over 
impact parameters gives 

J d2bx2„(s,&;p9 = <.(s;Pt), (21) 

where (T|^(s;pt) is the integrated cross section for pro- 
ducing exactly n-dijet pairs. [The X2n(s, 6; Pf) correspond 
to the "exclusive" cross sections of Ref . [20] , the quotes 
referring to the fact that these cross sections are still in- 
clusive in soft fragments.] 

Now we reconstruct the hard dijet contribution to the 
inelastic profile function by writing down the expression 
for the inclusive cross section for fc-dijet production in 
terms of the exclusive cross sections: 

X2k{s,b;pf)^Yl JX2n(5,&;K)- (22) 

n>k ^ ^ 

The combinatorial factor counts all the ways an n-dijet 
event can contribute to the inclusive /c-dijet cross section. 
Equation can be inverted to obtain the exclusive 
impact parameter cross sections in terms of the inclusive 
ones, 

X2fc(s,&;p9 = E (l)i-^r''x2nis,b;pt). (23) 

n>k ^ ^ 

The total inelastic profile function is obtained by sum- 
ming all the exclusive components. Using Eq. ([23l) . we 
obtain 

oo 

rf^Us,b;p^)=Y,X2k{s,b;pd^ 

oo 

= E(-i)""'^2„(s,&;p^). (24) 

n=l 



In principle, if all the inclusive n-dijet impact parameter 
dependent cross sections X2nis,b;p^) are known, then it 
is possible to obtain exactly the dijet contribution to the 
inelastic profile function by summing all the terms in the 
last line of Eq. p4|. In practice, higher orders in n need 
to be modeled or approximated. 

Consistency requires the hard dijet contribution to the 
total inelastic cross section to be less than the actual 
total inelastic cross section so 

r^Us,b;p'l)<TZfUs,b), (25) 

where the right side is the "actual" inelastic profile 
function, which could be obtained from either a mea- 
surement or a model extrapolation. Hence, Eq. (p5)) 
provides a means of checking that an expression for 
rj{[J^jg(s, 6;Pj), constructed from Eq. is consistent 

with other methods for obtaining the total inelastic pro- 
file function. A violation of Eq. (|25t means either that 
the model/extrapolation is incorrect or that there is a 
problem with the X2n(s, &;Pt) used in Eq. (j24p . 



IV. UNCORRELATED SCATTERING 

The simplest and most common way to obtain an ex- 
plicit expression for rj{|J^jg(s, 6;pj) from Eq. (|24)) is to 
assume that all partonic collisions occur completely in- 
dependently from one another. It was shown in Ref. (27j 
that the inclusive impact parameter dependent cross sec- 
tion for production of n dijets is then 

X2n{s,b;pt) = ^X2{s,b;pir- (26) 

This can be inserted into the last line of Eq. and 
summed to reproduce the familiar unitarized eikonal-like 
expression, 

oo 

r^,^{s,b;p^) = -{-ir-'x2{s,b;pir = 

1 ^* 

n—1 

= l-exp[-x2(s,6;pD] • (27) 

The single, double, and triple scattering terms are rep- 
resented graphically in Fig. [31 (This kind of graphical 
representation will be useful later for describing combi- 
natorial factors when correlations are included.) Each 
circle-cross represents a hard scattering event. The un- 
correlated assumption of Eq. (j26p is symbolized by the 
absence of any lines connecting the different hard colli- 
sions - each graph is simply X2(s, b) raised to the appro- 
priate power. 

For many practical purposes, Eq. (P7)) is sufficient. In 
general, the reconstructed profile function simply needs 
to reproduce the correct pQCD expression at large b 
where multiple collisions are very rare, while the minimal 
unitarity requirement that rjj"?ijg < 1 should be enforced 
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FIG. 3: Graphical representation of the term in the series for uncorrelated scattering - 
of Eq. (|27p . assuming no correlations. Spectator partons are not shown. 



the first three terms in the second line 



at small b. That is, the basic requirements in the high 
energy limit are, 



-pincl 



dijots('5: b) 

inel , 
dijets 



X2(s,6;Pt; 
1. 



(28) 
(29) 



In the high energy limit, rjf|j|ljj,(s, &) is expected to ap- 
proach one at small b (black disk limit). Equation (P7)) 
is completely satisfactory as far as conditions (|28l29p are 
concerned. For large 6, only the first term in the series in 
Eq. ()27|) — i.e. single scattering — is important. How- 
ever, in a precise treatment, one should also account for 
potential for violations of the consistency requirement 
in Eq. psp in the high energy limit and at intermedi- 
ate values of b where corrections of order X2('S, b;p1)^ are 
non-negligible. Indeed, such consistency problems were 
found in Ref. \M- 



V. ORGANIZING CORRELATIONS IN 
MULTIPLE COLLISIONS 

A. Impact Parameter Independence 

Deviations from uncorrelated scattering can arise from 
multiple sources. As discussed in the Sec. U correlations 
can be generated both in perturbative evolution equa- 
tions and in nonperturbative models. 

Correlations will also be induced by kinematical con- 
straints. We will assume, however, that most active par- 
tons have small enough x that these constraints are unim- 
portant, at least for the first few terms in the series in 
Eq. (I?fl) . For this paper, we will assume that the in- 
coming partons that take part in multiple hard collisions 
move nearly parallel with transverse momentum of order 
~ 1/Aqcd- That is, they have momentum typical for 
bound constituents of the incoming hadrons. In general, 
if is allowed to be larger than a few GeV, the partons 
will undergo DGLAP evolution, and hence may include 
partons with larger transverse momentum. Furthermore, 
one expects significant dependence of CTe// on the hard 
scale at large [22|. In such cases, it is possible that 
correlations may be understood as arising from parton 
evolution. However, conflicts with Eq. (|25p become less 
likely at larger p^. 



Therefore, we organize the description of correlations 
around the assumption that the effect is to introduce a 
simple (impact parameter independent) rescaling from 
the uncorrelated case. As a first example, we reconsider 
double hard collisions. Equation (pS)) gives the uncorre- 
lated expression 



X4(s,&;Pt) = ix2(s,&;Pt)^ 



(30) 



which should be replaced in the correlated case by 

X4(s,&;pD^ ^(1 + ^4(s))x2(s,6;p9', (31) 

where 774(5) parametrizes the deviation from uncorrelated 
scattering. Our strategy is to estimate the size of the 
double correlation by directly fitting Eq. (I5T|) to experi- 
mental data, given the constraint that X2(s, 6;Pt) is fixed 
by the GPD in Eq. (|19p . Note that we place no condi- 
tion on the b integral of Eq. (|30|) . In particular, we do 
not use the approximation in Eq. ([2]). In general, the 
correlation correction will also depend on both p^ and b. 
For our analysis, we will not explicitly write the pi argu- 
ments in Eq. pop because we are mainly concerned with 
correlation corrections in the limited range of pi where 
Eq. psp becomes problematic within the usual eikonal 
picture. As we will see, neglecting the b dependence in 
774(5) will allow for a direct parametrization of the cor- 
relation correction in terms of experimentally observed 
double scattering rates. It is likely that this is a very 
rough approximation, but it will allow for a first estimate 
of the role of correlations at large impact parameters. We 
also remark that the dynamics responsible for confine- 
ment are likely to induce large correlations regardless of 
impact parameter. We will discuss possible b dependence 
in greater detail in Sec. I VIII 

In experiments the effect of double partonic collisions 
is most commonly represented by the observable. 



2 ar{s;p1) ■ 



(32) 



In the uncorrelated case, using Eq. 
Eq. ^ yields 



1 



Jd^hP2{s,b;plY 



301) and Eq. (|T9| in 



(33) 
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FIG. 4: Graphical representation of the extra contribution 
(— l/2)774(s)x2(s, )^ due to double partonic correlations 
in Eq. (PT|). 



In general, the value of acs can be fitted to experimen- 
tally measured values by changing the width or shape of 
P2(s, b;p^). However, in our approach P2(s, b;Pf) is fixed 
by experimental measurements of the GPD, so the width 
of P2(s, ^;Pt) is not a free parameter. 

If Eq. pip is used in Eq. p2p , one obtains for (Toff in 
the correlated case 



(1 + 7^4(5)) /d2bP2(s,6;p^)2- ^ ' 

With P2(s, b'jPf) fixed by the two-gluon form factor, one 
can only tune Eq. ([34]) to the measured value of CTcff by 
adjusting 774(5). Note also that, although we have as- 
sumed impact parameter independent correlations, the 
relative rate of double collisions compared with single 
collisions, X2is,b)/x4:{s,b), still depends on impact pa- 
rameter. 



Eq. in]): 

[(1 + 774(5)) X2(s,&;p^)2] ^^^s,b;pl). (36) 

There is, therefore, an extra contribution equal to 
^^4^(^2(5, &))^ for each of the (2) = 3 ways a pair of 
incoming bound partons can become correlated. This is 
illustrated graphically in Fig. |SJ which shows the addi- 
tional contributions that must be added to the n — 3 
term in Eq. p7]) /Fig. l3l The expression for X6(s, b;p1) is 
therefore 

Xeis, b-p'i) = i (1 + 3r;4(s)) ^2(5, b-p'if. (37) 
6 

Following this example, it is now clear how to include 
double correlation corrections in n-parton scattering. In 
an n-parton collision, there are (2) additional contribu- 
tions equal to ■^'rii:{s)x2{s,b]p'i)'^ . In terms of diagrams 
like Fig.jS] this corresponds to all the ways that two colli- 
sions can be connected by a single zigzag line. Therefore, 
to include double correlations in the description of the in- 
clusive 7i-dijet cross section, the uncorrelated relation in 
Eq. ([26]) should be replaced with. 



Double Partonic Correlations in Multiple n\ \ L j 

Collisions (38) 



The effect of double correlations in an 77,-collision event 
may now be organized in a very convenient way. We 
start by looking at how r^"?^^g(s, &;pj) is modified by 
the inclusion of double correlations. In Eq. ([?!]) ■ 774(5) 
parametrizes the deviation of X4(sj^) from the uncorre- 
lated case, 774(5) — 0. It represents a correction to the 
assumption in Eq. ^ that the integrated double parton 
PDF is simply a product of the standard PDFs. The 
additional term proportional to 774(5) is represented by 
Fig.HJ The zigzag line connecting the two hard collisions 
may be thought of loosely as representing the effect of 
summing all soft gluons exchanged between the nearly 
parallel incoming and outgoing partons. We call 774(5) 
the double correlation correction factor. 

Next, we reconsider the uncorrelated description of 
triple parton scattering, graphically represented by the 
third term in Fig. [3] With no zigzag lines, we get the 
naive uncorrelated contribution from Eq. (|26p 

X^{s,b-pX)^\x2{sMv'tf- (35) 

For each pair of incoming partons there is another dou- 
ble correlation correction. In other words, for each pair 
of colliding partons there is another replacement like 



Using Eq. (|38]) in Eq. ([24]) and summing over all n pro- 
duces an analytic expression for the hard dijet contribu- 
tion to the total inelastic profile function, 

ri-'(5,6;p9 = 

E^^T^ U + ^?4(5) \ X2(5,6;p,T 
= 1-exp \-X2{sMvX)\- 
- ^ X2(5, b-pTi" exp [-X2(5, b-pT)\ . (39) 

Note that this equation respects the basic requirements 
of Eqs. (|28I29I) as long as 7^4(5) > 0. If 774 is allowed to be 
less than zero, then there is a potential for Eq. ([39| to be 
greater than one for some intermediate impact parame- 
ters. The third line is the standard unitarized expression 
Eq. ([??]) for the profile function, familiar from the eikonal 
model, while the last line is a correction due to double 
correlations. The double correlation correction term con- 
tributes a power of ^2(5, 6)^ in a series expansion in small 
X2(5,6;Pf). Therefore, it will only become important at 
impact parameters which are small enough that terms 
proportional to X2(5,&;Pt)^ are non-negligible. 
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FIG. 5: Graphical representation of the extra contribution to the n — 3 term of Eq. (|27|l due to double partonic correlations. 



C. Higher Correlations 

In the last section, we assumed that only a single pair 
of partons can become correlated. This is reasonable if 
the goal is simply to account for double correlations at 
large or intermediate impact parameters where the con- 
tribution from double collisions [order X2{s,b)'^] is sig- 
nificant, while contributions from triple collisions [order 
X2 (s, b)^] and higher are negligible. The explicit sum over 
all collisions in Eq. p9|) is needed to produce an analytic 
expression for the corrected inelastic profile function that 
still satisfies Eqs. (|28l29p and is less than unity for all b. 

In reality, there are of course corrections from triple 
and higher correlations; in our graphical representation, 
triple correlations are represented by zigzag lines con- 
necting three of the interaction points — see Fig. [6l The 
contribution from triple correlations becomes important 
only at order X2is,b;p'^)^. If powers of X2is, b;Pf)^ are 
significant, then we can iterate the steps of Sec. IV Bl by 
replacing Eq. (|37)) with 



1 + 3?74(s) + %(s))x2(s,6;p9^ (40) 



exactly analogous to Eq. ((3T|) for double correlations. 
The rjQ^s) parametrizes the correction from triple cor- 
relations in triple and higher partonic collisions. 

Taking into account the contribution from triple corre- 
lations to n > 3 collisions, and including the appropriate 
combinatorial factors by counting all ways of connecting 
three hard collisions, we then recover Eq. (15^ . but with 
a triple correlation correction term equal to 
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X2(s,&;Pt) exp [-X2(s,&;Pt)] • 



Now it is a simple matter to generalize the steps from 
Sec. lVBl to the arbitrary case of n correlation corrections. 
The resulting general expression for the inelastic profile 
function is 



= l-exp[-X2{s,b;pt)]- 
- E X2is, b;pir exp [-X2is, b;pf)] . 

n=2 

(41) 

The series after the first line includes all correlation cor- 
rections. There is a new correlation correction factor 



V2nis) for each number n of collisions. In the series repre- 
sentation for rj|i;jg(s, 6;Pt), a correlation correction factor 
r]2j{s) is accompanied by powers X2is, b]piy or higher. 

Equation (|¥T]) always satisfies the basic conditions 
Eqs. (|28I29|) of a unitarized profile function. [Although, 
depending on the signs of the ?72n(5), it may need to be 
checked that the profile function does not exceed unity 
for some intermediate value of b.] Contributions from 
higher n correlation corrections are suppressed by fac- 
tors of {x2{s,b)y''/n\ and can be neglected so long as 
X2{s,b) is sufficiently smaller than 1. By truncating the 
series at larger n, we obtain an increasingly refined de- 
scription of the b tail at moderate to large b. By using 
models of multiple collisions to obtain the ?72n(s), or di- 
rectly parametrizing the size of correlation corrections 
from experimental data, it should therefore be possible 
to reconstruct an inelastic profile function that respects 
Eq. dig). 

Unfortunately, there is as yet very little direct exper- 
imental knowledge of r]2nis) for n > 2. However, in the 
next section we will argue that even when only double 
correlations are included, the corrections are important 
at moderate to large impact parameters. Once data are 
available, steps analogous to those in Sec. IV Al can be 
used to parametrize rjQ^s). As in Ref. @, we can define 
the triple effective cross section. 



, T.2 ^ l aT'{s-ptf 
Then, including up to triple correlations, we have 



1 



(1 + 37?4(.S) + 776(.S)) / <Ph P2(S, b-plf 



(42) 



(43) 



From Eq. (j43p , we can calculate the correction from dou- 
ble correlations to cr^ with triple correlations neglected. 
In the next section, we will find values of 1.3 or 2.1 for 
774(5) at currently accessible energies and pi ~ 2.5 GeV. 
These give ct^ — 12.8 mb and 10.5 mb respectively, com- 
pared with — 28.3 mb for the case with no double 
correlations. 

We end this section by pointing out that, in princi- 
ple, the steps leading to Eqs. P9I41I) remain valid if we 
allow the r]2n{s) to have impact parameter dependence. 
The hypothesis of impact parameter independent corre- 
lations is only needed if we wish to estimate the size of 
correlations by using Eq. (IMl) . 



9 



FIG. 6: Graphical illustration of triple correlations in triple parton scattering. 



VI. NUMERICAL ESTIMATES 

A. Total Cross Sections 

Any numerical results that we obtain using Eq. (j41[) 
should be compared with common extrapolations of the 
total inelastic profile function to high energies so that 
consistency with Eq. ([25]) can be verified. A standard 
parametrization of the total profile function takes the 
form 

with B{s) ^ Bq + a' Ins. Regge theory fits give a ~ 
0.25 GeV-2 for the rate of growth of B{s). For the LHC 
energy of ^/s = 14 TeV, a survey of common models 
and extrapolations in the literature [29l - [33j suggests the 
following as a range of reasonable parameters: 

90mb < crtot(yi= 14 TeV) < 120mb (45) 
19GeV"^< B(VJ=14TeV)< 23GeV"^. (46) 

For example, in Ref. it is found that Eq. (gl]) with 
r(6 = 0) = 1 and B = 21.8 GeV'^ is in very close 
agreement with the Regge parametrization of Ref. [s^ as 
well as with the non-Gaussian model of Ref. [s^. [Some 
fits put the maximum from the total upper error band for 
(ytot{\fs = 14 TeV) at around 130 mb. However, a total 
cross section this large would also require a very large B 
to avoid having a profile function greater than unity at 
small b.] 

B. Correlated vs uncorrelated partons 

To calculate tTcff within the uncorrelated assumption, 
we use Eq. and obtain X2is,b;Pf) from the two- 

gluon form factor, as in Eq. (fT9|). We use irig « 1 GeV 
which works well for 0.03 < a; < 0.05 and large pt, 
relevant for most Tevatron data. More data on J/V' 
electro(photo)-production and deeply virtual Compton 
scattering in this range of kinematics would be very de- 
sirable for improving the accuracy of the determination of 
the b dependence of quark and gluon GPDs at a; ~ 10~^. 

The value of aes obtained from Eq. is then about 
34 mb. At small x, the width of the X2(s,&;Pt) grows 
and results in an even larger value for Ccff- The precise 



rate of growth of the radius at small x and fixed pt is not 
currently well established but will likely become clearer 
as new data become available. For the x dependence, we 
will use the parametrization in Ref. [iTj . 

The 34 mb calculation obtained with the uncorrelated 
approximation should be compared with the measured 
value of 14.5 mb [l| from the CDF collaboration taken 
at a center-of-mass energy of y/s = 1.8 GeV. The un- 
correlated calculation is roughly a factor of 2.3 too large, 
implying that it is unsafe to neglect corrections from cor- 
relations. At a minimum one should keep the 774(3) term 
in Eq. ([M)) with a correlation correction factor 774 « 1.3. 

It was argued in Ref. [ll] that the analysis in [l| ac- 
tually overestimates aes- If three-jet events are taken 
into account (to make the cross section truly inclusive), 
then a new estimate is CTeff ~ 11 mb. This suggests 
that the correlation correction is closer to 774 w 2.1. 
More recent measurements from the collaboration 
find (Toff = 15.1 mb [3], without cuts on three-jet events. 
So the precise size of a^s remains unclear. We remark 
that (Tcff may depend on pi, which may lead to differences 
in measured values 22] . 

The CDF measurements in Jl] find that correlations 
depend weakly on x, suggesting that 774 may be roughly 
constant with energy. Therefore, we will test the effect 
of using 1.3 < V-iiV^ = 14 TeV) < 2.1 in the calculation 
of the inelastic profile function using Eq. p9l) . Plots of 
Eq. ([5^ are shown as dotted curves in Figs. [7Ja-b). In 
these calculations we allow mg to vary slowly with x and 
P( in accordance with the parametrization in Ref. [l7| . 
The total inclusive dijet cross section is the same as what 
is used in fisl. b ased on the CTEQ6M gluon distribu- 
tion function [33] with a K factor of 1.5. With next-to- 
leading-order PDFs being used, the K factor is closer to 
1.2. However the dominant contribution to final states 
typically involves at least three jets, corresponding to 
K — 1.5 for our calculation [s^. See also the discussion 
in Ref. [18]. 

In Fig.lTJa), we have used 7/4 1.3 while in Fig.[71[b)we 
have used 7/4 w 2.1. In addition, we have tested the 
sensitivity to higher correlations by including terms up 
to n = 4 in Eq. (HJ), and using the approximation 778 « 
776 774 = 77. The resulting curves are shown as dashed 
lines in Figs. [7Ja-b). The suppression at large b from 
double correlations could in principle be spoiled if the 
triple correlation is large and positive, so we have checked 
the case where r/Q > 0. Then, to avoid the possibility that 
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r > 1 at very small b, 773 is made positive. 

For comparison, the completely uncorrelated eikonal- 
type expression, Eq. ([ST]) , is plotted as a dash-dotted 
curve in Figs. [71[a-b). Note that there is a substantial 
difference between the correlated curves and the uncor- 
related expression for 0.8 fm< 6 < 2 fm in both plots. 
In all cases, the correlations result in a suppression of 
the total inelastic cross section from dijets [calculated by 
integrating r^"j°j';jg(s, 6) over b] by more than 15%. The 
shaded regions mark the area covered by the standard ex- 
trapolations of the total inelastic profile function. They 
correspond to Eq. (j44p with the range of parameters in 
Eqs. (|45l |46)) . The pt cutoff in all cases is fixed at the 
typical value of « 2.5 GeV. 

The uncorrelated curve lies entirely above the shaded 
area for & < 1.6 fm, in violation of Eq. That is, 

the hard contribution to the total inelastic cross section 
is larger than the total inelastic cross section itself for 
much of the essential range of impact parameters. The 
curves that include double or quadratic correlations ex- 
hibit greater consistency for the full range of b for both 
Figs. [Tfa) andjTfb). In the case of the moderate sized 
correlation corrections in Fig.lTJa), the effect of including 
triple and higher correlation corrections is rather small 
compared with the case where only corrections from dou- 
ble correlations are kept. Including higher correlation 
corrections does seem to smooth out the shape of the 
profile function. (We have also assumed all correlation 
corrections to be positive.) However, the higher correc- 
tion terms only become significant at small impact pa- 
rameters where the profile function is already close to 
unity anyway. If the correlation corrections are larger, as 
in Fig.[7l^b), then the higher n> 2 correlations are more 
significant. 

Now let us consider what is needed for the radius of 
the hard overlap function if all the r?2n(s) are set to zero 
in Eq. (PT|) (reducing to the standard eikonal formula). If 
a small value of is used to evaluate apQQj^{s;p'j:), then 
fits of the total cross section to current data require a very 
narrow width for the overlap function [l2l|. In theoret- 
ical calculations, a narrow overlap function is obtained, 
for example, in the semiperturbative approach proposed 
in [3I] where the radius of the hard overlap function de- 
creases with energy. In Pythia the hard overlap function 
is modeled by the double Gaussian parametrization Q 



-62 



(1 - 0)^ 
^^(^)-S^^^Pl2«- 

2/3(1-/?) 
H 2 , 2 6xp 



^expj^^ (47) 



with 02 = 0.4 fli and /? = 0.5 (in Tune A). The radius in 
Eq. ()47p does not vary with energy. We determine ai by 
using Eq. (|T7|) in Eq. ([55]) for acS with uncorrelated mul- 
tiple hard scattering and fixing it to the measured CDF 
value. In Fig. [9] we show Eq. (|47l) with ai calculated us- 
ing (Jeff = 14.5 mb (Ref. 1]) and doff = H mb (Ref. 
For comparison we have also plotted the overlap func- 



tion Eq. (1161) obtained from the two-gluon form factor at 
y/s = 14 TeV, and the overlap function from Ref. [sl ]. 
also at y/s = 14 TeV [1^. A comparison of these curves 
illustrates that models of the overlap function tend to be 
much narrower than what is expected from the two-gluon 
form factor. 

Because ai and 772 are fixed by the measured value of 
CTcft, Eq- (|47l) yields by construction the same rate of dou- 
ble collisions as Eq. (jSH). For triple and quadruple colli- 
sions the two approaches give roughly similar rates (same 
orders of magnitude). Therefore, in practical situations, 
using Eq. (|47p with a narrow peak may be an economical 
way to model the effects of correlations. However, our 
basic aim in this paper is to incorporate the maximum 
amount of available experimental input into the descrip- 
tion of hard collisions by using the factorization theorem 
and parametrizations of the GPD to describe the overlap 
function. Using the overlap function obtained directly 
from the t dependence of the J/i/" photoproduction cross 
section requires either that a larger transverse momen- 
tum cutoff {p1 > 3.5 GeV) is used, or that double corre- 
lations are incorporated by using Eq. ([39| with ri2n > 0. 
Otherwise, there are potential problems with the consis- 
tency relation Eq. (pS)) . even for a relatively large b. 

One source of uncertainty is the shape of the overlap 
function P2{s, b]p1). A Gaussian form, for instance, may 
be preferred to Eq. (|16p . Therefore, we have repeated 
the calculation of Fig. [71 but now with 



P2{s,b;p't) 



1 



27rfeo 



exp 



26n 



(48) 



The parameter 60 is fixed by requiring that the average 
b\ 



(62) = / d^hb^P2{s,b;pt) 



(49) 



is the same for both Eq. (|48p and Eq. (1T51) . The resulting 
plots are shown in Fig. [51 The drop with b is slightly 
steeper at intermediate b in Fig. [3 but otherwise the 
plots are very similar. We also point out that a recent 
experimental study ;39] finds good agreement between p 
and electroproduction data and the dipole form for the 
two-gluon form factor. 

Another source of uncertainty is the contribution from 
diffraction to the inelastic cross section, which is expected 
to be much more peripheral than generic inelastic interac- 
tions. This is known already from analyses of the diffr ac- 
tive processes at lower energies (ioj and should be even 
more prominent at ^/s = 2 TeV and above where inelas- 
tic diffraction cannot occur at small impact parameters, 
and where the interaction is practically black. Inelastic 
diffraction constitutes a significant fraction of the inelas- 
tic cross section at -^i = 2 TeV, 25% - 30 %, and is 
expected to remain significant at the LHC. 

Hence, in the region where we use the consistency re- 
quirement, Eq. (I?51) . a large fraction of r'"°' is due to 
inelastic diffraction. At the same time the Tevatron data 
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FIG. 7: Inelastic profile functions calculated with and without correlations for ^/s = 14 TeV and pt = 2.5 GeV. The shaded 
region corresponds to the range of typical extrapolations. The dash-dotted curve corresponds to the standard eikonal expression. 
The dotted curve is the inelastic profile function including the double correlation correction in Eq. (|39} with (a) 774 = 1.3 [ll] 
and (b) 774 = 2.1 [s^. The dashed curve is with the triple quadruple and correlation corrections from Eq. (|4ip using r) — tia — 
V6 = Vs =(a)1.3 and (b)2.1. 



on the jet production in diffraction indicate a large sup- 
pression of jet production as compared to expectations 
based on the use of the diffractive PDFs measured at 
HERA. This indicates that for the values of which 
we discuss the fraction of the inelastic diffractive events 
with jets is significantly smaller than 1. Hence, Eq. ([^5)) 
should be applied to r''''\b) - V^'^^b). This implies that 



an even larger pt cutoff is needed. A more quantitative 
analysis of this effect requires modeling of the inelastic 
diffraction profile function and of the dynamics of the 
suppression of jet production in diffraction. We leave 
such investigations to future work. 
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FIG. 8: Inelastic profile functions calculated exactly as in Fig. |S] but now with the Gaussian form for the overlap function. 



VII. TESTS OF IMPACT PARAMETER 
DEPENDENCE 



A large source of uncertainty is in the role of b depen- 
dence in the correlation corrections. As far as we know, 
there are currently no predictions of the impact param- 
eter dependence of nonperturbative correlations in mul- 
tiparton distributions. If all correlations are localized at 
small impact parameters, then at large impact param- 
eters one simply recovers the uncorrelated model. We 
find such scenarios unlikely, however, since the binding 



interaction between any constituent partons should be 
expected to be large, regardless of impact parameter. By 
relaxing the assumption of b independence for all corre- 
lation corrections, it is possible to reconstruct arbitrarily 
different shapes for the profile function, though in princi- 
ple this arbitrariness can be reduced by future measure- 
ments of higher correlations. In Figs.lTJb) andjSJb), the 
dip at intermediate b that occurs when only double cor- 
relations are included suggests that higher correlations 
should be included. However, a smooth form for the in- 
elastic profile function can also be recovered if we allow 
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FIG. 9: Overlap functions obtained from the two gluon form factor, Eq. (|16p (solid curve), the Godbole-Grau-Pancheri- 
Srivastava (GGPS) model Ref. [s^ (dashed curve), and the Pythia Tune A overlap function Eq. (|47p fitted to a^fi — 14.5 mb 
(dash-dotted curve) and a^a = 11 mb (dotted curve). 



for modest impact parameter dependence in the double 
correlation correction. As an example, we use instead of 

774 = 2.1, 

mib) - e-'"' + (50) 

This gives numerically the same cTeff as Eq. p4p with 
774 = 2.1, but now with correlations peaked at small b and 
with a weakly falling tail at large b. The plot analogous to 
Fig. [TJb) is shown in Fig. [TUl The curve with only double 
correlations is now closer to the corresponding curve in 
Fig. El^b) which was for a larger a^f / . However, the main 
conclusion of the previous section remains valid — that 
correlation corrections of size roughly 1 to 2 are needed 
even at relatively large b in order to have consistency 
with Eq. (|25l) . 

It is possible to visualize why an enhancement in cor- 
relations at large impact parameters is likely by consid- 
ering the following simple model: Consider scattering at 
a large impact parameter |b| = p and assume that col- 
lisions are between the pion clouds in each nucleon with 
the core of the other nucleon. Then the cross section for 
scattering in a range of impact parameters from p — r^r 
to p + r^r (where is the pion radius) is proportional to 
the cross-sectional area of the pion and the probability 
prob^(/9) to scatter if the pion cloud from one nucleon 
overlaps with the nucleon core of the other: 

dfJcor oc TTrl X prob^(p). (51) 

By contrast, without correlations the cross section is pro- 
portional to the area of the annulus between p — r^^ and 
/o -|- and a different probability, 

rfcTuncorr OC 47rpr^ X Y>Yoh^^^^^^{p) . (52) 



Requiring that these two expressions give the same cross 
section means that 

The cross sections for double scattering are the areas 
times of the square of the probabilities for single scat- 
tering. The ratio of the double scattering cross sections 
obtained for the correlated and uncorrelated cases is then 

7rr^[prob^(p)]^ ^ 4p ^^^^ 
4p7rr^[prob^,„,„„(p)]2 r^' 

Hence, in this simple picture one expects the true double 
scattering cross section to be enhanced at large p relative 
to what is expected if correlations are neglected. 

A natural question is how to test the dependence of 
correlations between partons as a function of b. One pos- 
sibility is to study the dependence of multijet production 
on the associated hadron multiplicity away from the ra- 
pidities and angles where hadron production due to the 
fragmentation of jets is important. The distribution in 
b of events with dijets is given by the overlap function 
Piib^ [Eq. [TH]. Let us now consider the distribution over 
the accompanying multiplicity Pm{N/{N)) where N is 
the observed hadron multiplicity, and (N) is the average 
multiplicity in minimum bias nondiffractive events (the 
CDF collaboration reported an average multiplicity for 
events with a Z boson trigger for angles where gluon ra- 
diation effects associated with the production of Z are 
small by a factor ^ 2 larger than in the minimal bias 
inelastic events j4l|). 

Both soft and hard interaction mechanisms of hadron 
production lead to a monotonic increase in the average 



14 



0.6 




FIG. 10: Curves analogous to Fig. [7jb), but now using the 6-dependent correlation correction in Eq. (fSOl 



accompanying multiplicity with decreasing b. In the first 
case this is due to an increase of the contribution of the 
multi-Pomeron exchanges; in the second it is due to an 
increase in the probability of multiple (semi)hard inter- 
actions. Hence, to a large extent fluctuations in the mul- 
tiplicity are due to the distribution of collisions in b. To 
simplify the following discussion we will neglect fluctua- 
tions of multiplicity at a given b, though one can model 
these effects as well using current Monte Carlo models. It 
seems natural to expect that such fluctuations will smear 
resolution in b but will not remove the generic trend of 
decrease of the multiplicity with increase of b. 

Under these assumptions, we can identify intervals of b 
that correspond to production of i hadrons not belonging 
to the hard collision (this could be, for example, hadrons 
at the central rapidity interval not affected by hadron 
production in the hard process leading to production of 
dijets). 



P2{b)d%^p, 



(55) 



where pi is the probability for producing exactly i 
hadrons. Given values oipi and P2ib), we may construct 
the number of hadrons N{b) corresponding to a given 
impact parameter: 



N{b) 
Nib) 



iorbe [b2,bi] 
for 6 € [6.3, &2] 



(56) 
(57) 



One can write similar equations for the higher moments 
of N^^\ though in this case sensitivity to the fluctuations 
at a given impact parameter becomes larger. 

If correlation corrections are impact parameter inde- 
pendent, then P^ib) is simply proportional to the square 
of P2{b)- However, if correlations are concentrated at 
small impact parameters, then Pilb) is more sharply 
peaked at small b and falls off faster at large b. Then the 
integrand in Eq. (|58)) is more localized at small b where 
N{b) is large, and the average accompanying multiplic- 
ity for four-jet events will be larger than in the uncorre- 
lated case. A similar analysis extends to higher moments 
(A^^*^). One should emphasize here that the multiplicity 
of hadron production in the rapidity interval between 4 
jets may be affected by various effects of color correla- 
tions, etc.; hence it is desirable to look for the change of 
the multiplicity at the rapidities sufficiently remote from 
the region of 4 -jet activity. This is feasible for the LHC 
detectors with a good acceptance in a large rapidity in- 
terval. 

In addition, if the correlations are present at all im- 
pact parameters they should be manifested in the hard 
diffractive processes which correspond to scattering at 
large impact parameters. One could consider both cases 
of single and double diffraction with production of two 
and four jets: 

pp -^p + X{2jets + Y,Ajets + Y); (59) 
pp ~^pp + X{2jets + Y,Ajets + Y). (60) 



Next we can consider production of four jets. We may 
calculate the multiplicity of events with at least four 
jets by integrating the number of collisions N{b) over b, 
weighted by the probability density Pi{b) for a four-jet 
event: 



(Ar(4)) ^ / p^(h)N{b)cPb. 



(58) 



Correlations between the partons should also enhance 
the cross section of the exclusive channel when the light- 
cone fraction carried by two of the interacting partons of 
one of the nucleons is close to maximal: {xi-\-X2) / xp ~ 1. 
Such a contribution should be enhanced if —t is large 
enough (few GeV^) to squeeze the transverse size of the 
exchanged ladder (see Fig. [TT|) . 
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FIG. 11: Double Pomeron process with production of two 
pairs of dijets. 

VIII. DISCUSSION AND CONCLUSIONS 

The main conclusion of this paper is that, given the dis- 
tribution of hard partons known from the gluon GPD, the 
hard contribution to the total inelastic profile function 
should probably be modeled using Eq. or Eq. (|1T|) 
with 1 ^ 77 ^ 2 rather than the usual one-minus- 
exponential shape that arises in a purely eikonal treat- 
ment. This avoids conflicts with general expectations for 
the total inelastic profile function (essentially a unitarity 
problem), when a small or fixed pi is used for the inclu- 
sive dijet cross section. The correlation corrections stem 
from a breakdown in the factorization ansatz of Eq. 

There are many uncertainties associated with the size 
of correlations at various impact parameters. Neverthe- 
less, we believe there is ample evidence that large corre- 
lations are important even at relatively large impact pa- 
rameters. Our sample calculations illustrate how includ- 
ing correlation corrections of roughly the size suggested 
by available experimental data can lead to greater con- 
sistency with the total inelastic pp cross section. Hence, 
questions of unitarization and consistency at small pi 
should be organized around a more precise determination 
of higher correlation corrections, with either experimen- 
tal input or theoretical modeling. 

As more data become available, and the range of al- 
lowed parameters narrows, it will be possible to use 
Eq. (1411) to obtain an increasingly refined picture of the 
profile function. More studies of the x dependence and 
rate of growth of the hard overlap function are needed. 
Furthermore, it will be important to establish, through 



measurements or theoretical models, the nf and energy 
dependence of the correlation corrections [221 ■ Although 
we have assumed that correlation corrections are impact 
parameter independent in this paper, it is possible to 
incorporate any impact parameter dependence into the 
same basic framework simply by allowing the r]2n{s) to 
depend on impact parameter. We have suggested pos- 
sible ways of testing for impact parameter dependence 
of correlations in Sec. I VIII It would also be very useful 
to have direct measurements of contributions from triple 
and higher correlations. For a recent discussion of this 
possibility for triple correlations at the LHC, see Ref. 
Furthermore, it was recently argued in Ref. [42] that dif- 
ferent numbers of collisions contribute incoherently when 
they are identified by their topologies. 

Fits to the total cross section can be obtained within 
the more common eikonal approach by using a very nar- 
row hard overlap function. However, the contribution to 
the total cross section arising from hi gh pt J ets depends 
on the blackness of the hard collisions [13, [3J| , and hence 
on the width of the hard overlap function. Therefore, it is 
phcnomenologically important to use the correct radius 
in the description of hard multiple collisions. In addition 
to being necessary for the construction of realistic simu- 
lations, sorting out these issues has the potential to lead 
to an improved understanding of the transverse structure 
of the proton at high energies. 
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